Genome-wide identification of DNA methylation provides insights into the association of gene expression in rice exposed to pesticide atrazine

Atrazine (ATR) is a pesticide widely used for controlling weeds for crop production. Crop contamination with ATR negatively affects crop growth and development. This study presents the first genome-wide single-base-resolution maps of DNA methylation in ATR-exposed rice. Widespread differences were identified in CG and non-CG methylation marks between the ATR-exposed and ATR-free (control) rice. Most of DNA methyltransferases, histone methyltransferases and DNA demethylase were differentially regulated by ATR. We found more genes hypermethylated than those hypomethylated in the regions of upstream, genebody and downstream under ATR exposure. A stringent group of 674 genes (p < 0.05, two-fold change) with a strong preference of differential expression in ATR-exposed rice was identified. Some of the genes were identified in a subset of loss of function mutants defective in DNA methylation/demethylation. Provision of 5-azacytidine (AZA, inhibitor of DNA methylation) promoted the rice growth and reduced ATR content. By UPLC/Q-TOF-MS/MS, 8 degraded products and 9 conjugates of ATR in AZA-treated rice were characterized. Two of them has been newly identified in this study. Our data show that ATR-induced changes in DNA methylation marks are possibly involved in an epigenetic mechanism associated with activation of specific genes responsible for ATR degradation and detoxification.